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After presenting a brief review of how simulations of QCD with imaginary chemical 
potential can be used to extract physical results, we analyse the phase structure of 
QCD with four flavours of dynamical fermions in the finite temperature - imaginary 
chemical potential plane, and discuss perspectives for realistic calculations. 



1 Introduction 

/ ^QCD \ 

The zero density QCD partition function, Z{V,T) ~ Trie r with 

Hqcb the QCD Hamiltonian, can be discretized on an euchdean lattice with 
finite temporal extent t = 1/T 

Z = J (l?C/I?7/>PV5)e-^^'=[^le~^'^[^^'^^'^l =y (I?C/)e-^^<=[^ldetM[C/] (1) 

where U are the gauge link variables and tp and ^ the fermionic variables, Sg 
is the pure gauge action, Sp is the fermionic action which can be expressed 
as a quadratic form in the fermionic fields in terms of the fermionic matrix 
M[U], Sf = i^M[U]ij. 

To describe QCD at finite density the grand canonical partition function, 

Z{V,T,fi) = Trie t j, where N = J c?'^a;i/'V is the quark number 
operator, can bcjjsed. The correct way to introduce a finite chemical potential 
fi on the lattice a is to modify the temporal links appearing in the integrand 
in Eq. as follows: 

Ut e°'^Ut (forward temporal link) 

Ul e^^^U} (backward temporal hnk) , (2) 

where a is the lattice spacing. Whilst Sg is left invariant by this transfor- 
mation, detAf[?7] is not and gets a complex phase which makes importance 
sampling, and therefore standard lattice MonteCarlo simulations, unfeasible. 
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The situation is different when the chemical potential is purely imaginary. 
This is implemented on the lattice as described in Eq. but in this case 
[/j e'^t^'Ut, Uj e-"'^"U}. This is like adding a constant U{1) back- 
ground field to the original theory; det M[U] is again real and positive and 
simulations are as easy as at /i = 0. 

The question then arises how simulations at imaginary chemical potential 
may be of any help to get physical interesting information. 

One possibility is analytic continuation □. Z{V, T, fi) is expected to be an 
analytical even function of /i away from phase transitions. For small enough 
fi one can write: 

log Z{^) ^ ao + a2/i^ + aifi'^ + 0{fi^) (3) 
log Z{fii) = ao- 02/1/ + a4M/ + 0{fi^j) . (4) 

Simulations at small /i/ will thus allow a determination of the expansion 
coefficients for the free energy and, analogously, for other physical quanti- 
ties, which pan be cross-checked with those obtained by standard reweighting 
techniques oQ. This method is expected to be useful in the high tempera- 
ture regime, where the first coefiicients should be sensibly different from zero; 
moreover the region of interest for present experiments (RHIC, LHC) is that 
of high temperatures and small chemical potential, with ^/T OA. This 
method has been already investigated in the strong coupling regime □, in the 
dimensionally reduced S-d QCD theoryO, and in full QCD with two flavoursB. 
The Taylor expansion coefficients can also be measured as derivatives with re- 
spect to /i at /i = 13'EI. 

Z(V,T, j/i/) can also be used to reconstruct the canonical partition func- 
tion Z{V, T, n) at fixed quark number n u, i.e. at fixed density: 

ZiV,T,n) = Tr ^ie-^S{N - n)) = ^Tr ^e'^ fj AOe^'^^^^) 



27r 



Aee-'"'^Z{V,T,ieT) . (5) 



1 

As n grows, the factor e~*^" oscillates more and more rapidly and the error in 
the numerical integration grows exponentially with n: this makes the applica- 
tion of the method difficult especially at low temperatures where Z{V, T, ifxi) 
depends_very weakly on /i/. The method has been applied in the 2-d Hubbard 
model E£l, where ZiV^ T, n) has been reconstructed up to n = 6. 

The phase structure of QCD in the T ~ ifij plane is also interesting by 

e'^^e —], Z{9) is 

clearly periodic in 6 with period 27r and a period 27r/3 is expected in the 
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confined phase, where only-physical states with N multiple of 3 are present. 
However it has been shown El that Z{9) has always perkjd 27r/3 for any physical 
temperature. Moreover the suggestion has been made a, based on a calculation 
in the weak coupling approximation, that discontinuities in the first derivatives 
of the free energy at = 27r/3(fc + 1/2) should appear in the high temperature 
phase. This suggests a very interesting scenario for the phase diagram of QCD 
in the T - ifj,i plane which needs confirmation by lattice calculations. 

We have recently started a program of simulations of QCD at finite imag- 
inary chemical potential, with both the aim of determining the phase diagram 
of QCD in the entire T - i^j plane, and exploring by analytic continuation 
the high T - small real chemical potential region. 

We have studied QCD with four degenerate staggered flavours of bare 
mass a - m = 0.05 on a 16** x 4 lattice, where the phase transition is expected 
at a critical coupling /3c — 5.04 Ell (the two flavour case has been studied 
as well eI). The algorithm used is the standard HMC $ algorithm. We will 
present here onlyA^ubset of our results and analysis. A complete presentation 
will appear soon t3. 

2 Results 

In order to understand the phase structure of the theory, it is very useful to 
look at the phase of the trace of the Polyakov loop, P{x). Let us parametrize 
P{x) = |P(S')|e*'^, and let (</>) be the average value of the phase. In the pure 
gauge theory the average Polyakov loop is non zero only in the dcconflned 
phase, where the center symmetry is spontaneously broken and (0) — 2fc7r/3, 
k ~ —1,0,1, i.e. the Polyakov loop effective potential is flat in the conflned 
phase and develops three degenerate minima above the critical temperature. 
In presence of dynamical fermions P{x) enters explicitely the fermionic deter- 
minant and Z3 is broken: the effect of the determinant is therefore like that 
of an external magnetic fleld which aligns the Polyakov loop along (0) = 0. 
In the high temperature phase the degeneracy is lifted: {(j)) = is the true 
vacuum and (</>) = ±27r/3 are now metastable minima. 

When Hi ^ 0, what enters the fermionic determinant is P{x)e^^, 6 = 
Hi/T, instead of P{x). Therefore the determinant now tends to align (</>) + 9 
along zero: its effect is like that of an external magnetic field pointing in 
the 6 direction. Hence one expects (</>) = —0 at low temperatures; at high 
temperatures the external magnetic field will still lift the Zz degeneracy, but 
now which is the true vacuum will depend on the value of 9. In particular 
one expects that for (fc — 1/2) < -^9 < [k -\- 1/2) the true vacuum is the one 
with (0) ~ 2fc7r/3 and that 9 = 2{k + l/2)7r/3 corresponds to first order phase 
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Figure 1 . Average value of the Polyakov loop phase as a function of the imaginary chemical 
potential for different values of /3. The vertical dashed line corresponds to 6 = fJ-i/T = n/3. 



transitions from one sector to the other: this is indeed the prediction of 
Roberge and Weiss. 

In Fig. we report our results for (</>) versus the imaginary chemical 
potential for different values of /3. Since T = l/{Nta) and A^f = 4 in our 
case, we have 9 ~ 4a/i/. For f3 — 4.94 and 5.01, which are below the critical 
P at fij — 0, /3c(m/ = 0) = /3c — 5.04, one has (0) ~ — = —Aafij, i.e. (0) 
is driven continously by the fermionic determinant. For f3 = 5.10, which is 
well above /3c we see that (0) ~ 0, almost independently of /i/, as long as 
6 < 7r/3, while for 6 > tt/3 there is a sudden change to {(/)) ~ —tt/3: we 
are clearly crossing the Roberge- Weiss phase transition from one Z3 sector 
to the other. At intermediate values, (3 = 5.065 and 5.085, {(/)) ~ until a 
critical value of a/^/, where it starts moving almost linearly with /i/ crossing 
continuously the Z3 boundary: in this case there is no Roberge- Weiss phase 
transition, but there is anyway a critical value of jij after which (0) is no more 
constrained to be ~ and can be driven again hy 9: as we will soon clarify, 
this critical value of ^/ corresponds to the crossing of the chiral critical line, 
i.e. the continuation in the T-fij plane of the chiral phase transition. 

We display our results for the chiral condensate in Fig. |^. We expect a 
periodicity with period 27r/3 in terms of 9. Moreover {ipip), like the partition 
function, is an even function of /i/: this, combined with the periodicity, leads 
to symmetry around all points 9 — rnr/S. with n an integer number, for {^ip) 
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Figure 2. Average value of the chiral condensate as a function of the imaginary chemical 
potential for different values of /3. The vertical dashed lines correspond to 8 = fii /T = 
{2k + l)7r/3. The continuous line in the lower picture is the result of a quadratic fit at small 
values of afij obtained at /3 = 5.10. 

as well as for the partition function itself. For f3 < (3c, (V'V') has a continuous 
dependence on afii with the expected periodicity and symmetries. For (3 > f3c 
the correct periodicity and symmetries are still observed but the dependence 
is less trivial. At /3 = 5.065 there is a critical value a/j,/ ~ 0.17 for which 
the theory has a transition to a spontaneously broken chiral symmetry phase: 
we are clearly going through the chiral critical line. The same happens for 
fi = 5.085 at am ~ 0.22: in this case we have proceeded further, observing 
also the transition back to a chirally restored phase at a/z/ ~ 0.30, which is, 
correctly, the symmetric point with respect to = 7r/3^. At /3 = 5.10 we never 
cross, when moving in /^j, the chiral critical line, but only the Roberge- Weiss 
critical lines ^. 

We have also performed runs at fixed /i/ and variable (3 to look for other 
locations of the chiral line iii_the T-fij plane, a detailed summary of all results 
will be presented elsewhere Il3. We present, in Fig. ^, a sketch of the phase dia- 



" We notice a discrepancy in the results obtained at the same /3 by multirewoighting tech- 
niques where the same symmetry cannot be observed. 

' Error bars for the determinations at /3 = 5.10 and on the critical lines (9 = n/3 and 
9 = it) arc probably underestimated. 
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Figure 3. Sketch of the phase diagram in the fii-fB plane. The filled cirles represents direct 
determinations of the chiral critical line location from our simulations. The rest of the chiral 
line has been obtained by interpolation and by exploiting the symmetries of the partition 
function. 

gram in the P-fii plane, as emerges from our data and by exploiting the above 
mentioned symmetries. We can distinguish a region where chiral symmetry 
is spontaneously broken (indicated as IV in the figure) and three regions (I, 
II and III), which correspond to different Z3 sectors and repeat periodically, 
where chiral symmetry is restored. The chiral critical line separates region 
IV from other regions, while the Roberge- Weiss critical lines separate regions 
I, II and III among themselves. In QCD with 4 staggered flavours, am — 0.05 
and fii — , the phase transition is expected to be first order EJ: assuming 
it continues to be first order also at fij ^ 0, we expect all the regions to be 
separeted by first order critical lines. 

It is interesting to illustrate the determination of the endpoint of the 
Roberg- Weiss critical line, (3e = 5.097(2). We have performed a simulation 
at exactly 9 = tt/S: starting thermalization at an high value of (3 from a 
zero field configuration we always stay on the border of region I, since on the 
16'^ X 4 lattice that we have used it is practically impossible to flip into region 
II through the critical line. As we decrease /3, always staying on the border of 
region I, we will meet the chiral critical line at (3 — (3e- Various quantities can 
obiouvsly be considered to signal this transition, but it is interesting to notice 
that in this case the baryon density can be taken as an exact order parameter. 
Indeed the baryon density, {b) = y-S-\nZ, is an odd function of /i, being 
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Figure 4. Imaginary part of the barion density as a function of fij for different values of /3 
(left— hand side), and as a function of /3 at 9 = fij/T = 2^ (right-hand side). 

Z an even function. Therefore, for an imaginary chemical potential, (b) is 
also purely imaginary and an odd function of fii. This, combined with the 
periodicity in /z/, leads to the expectation that {b){d = ^) — —{b){9 = 
The last relation clearly implies that (6) = at 6* = 7r/3, unless (&) is not 



continuous on that point. Thus a non-zero value of (6) = at 



implies 



the presence of the Roberge- Weiss critical line. On the right hand side of 
Fig. H the imaginary part of (6) at = ^ is plotted as a function of /?: one 
can clearly see a transition from a zero to a non-zero expectation value, which 
permits the determination of We have verified that the transition through 
/3_E is also visible in the chiral condensate: this implies that the Roberge- Weiss 
critical line ends on the chiral critical line. On the left hand side of Fig. ^ 
we present instead the imaginary part of (6) as a function of /i/ for different 
values of /? < /J^: in this case (6) is always zero and continuous at = ^, but 
it is interesting to note how it starts developing the discontinuity as /3 — > 

In order to translate results for /3c (/i/) into results for the physical crit- 
ical temperature, Tdp-i), we need the lattice spacing, a — a((3), in physical 
units. For instance using the values a(5.04) — 0.30(2) fm and a(5.097) = 
0.272(10) fm B we obain = 164(10)MeV and Te = 181(7)MeV. 

Finally we notice that for high temperatures and away from the critical 
lines, physical quantities show a clear non-zero dependence on the imaginary 
chemical potential, which is encouraging in starting the program of fitting 
the first terms of their Taylor expansion in fij and performing the analytic 
continuation to real chemical potential. As an example we have reported, in 
Fig. the result of a quadratic fit in /// for the chiral condensate a.t (3 — 5.10. 
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3 Summary and discussion 

We have clarified the phase structure in the imaginary chemical potential - 
temperature plane for full QCD with four staggered flavours. In particular 
we have confirmed the existence of the Roberge Weiss critical lines, located 
their endpoints and assessed their interplay with the chiral critical lines. 

We have checked that data at high temperature and small imaginary 
chemical potential can be safely fitted by a polinomial, when away from phase 
transition lines, making the analytic continuation to real /x feasible. 

The physical interesting region of high temperature and small chemical 
potential can now be studied by essentially three different techniques: 1) 
Direct calculations of derivatives; 2) Reweighting; 3) Analytic continuation 
from imaginary /x. Each method has its own merits and limitations, and cross 
checks among the three approaches should produce reliable results. 
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